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COLLIDING ELECTROMAGNETIC WAVES: 
DIAGONAL SOLUTIONS 



This chapter is a continuation of the review of all presently known exact 
solutions which describe the collision of plane electromagnetic waves, or a 
combination of both gravitational and electromagnetic waves. Attention 
is concentrated here only on diagonal solutions. These solutions may be 
considered as a generalization of the solutions representing the collision 
of gravitational waves with colinear aligned polarization that have been 
described in Chapter 10. It may be mentioned that the Bell-Szekeres 
solution that has previously been described in Chapter 15 also belongs to 
the class that is now being surveyed. 



18.1 The generation technique of Panov 

Panov (1978, 1979a) has obtained a class of exact solutions for colliding 
gravitational-electromagnetic waves using a generating technique due to 
Enss (1967). His method starts with a vacuum solution with aligned 
polarization, and can conveniently be described as follows. 

Taking the metric in the form of the Szekeres line-element (6.20), if 
a solution of (6.22) describing a collision of plane colinear gravitational 
waves is given by U Q , V Q , and M Q with W = 0, then a new solution with 
non-zero electromagnetic components is given by 



U = U Q , W = W = 

2 



= e Vo (cos 2 a + sin 2 a e ~^ Uo+v °^ 

o+v )y 



e -M = e -M a / 2 a + sin 2 ae -(U, 



(18.1) 



where a is an arbitrary constant. The Maxwell components are given by 

sin a cos ae~( u ° +v °^ 2 , rT Tr s 

$ = -*7 2 ^.2 _m +V) \ ( U °v + V ov) 

(cos 2 a + sin ae ( u °+ v °)) 

■ -(U+V),2 ( 18 - 2 ) 
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(cos 2 a + sin ae \ u °+ v °>) 
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It may immediately be seen from these equations that the boundary 
conditions for colliding plane waves are automatically satisfied for the new 
solution if they were satisfied for the original vacuum solution. The new 
solution therefore describes the collision of combined gravitational and 
electromagnetic waves with aligned polarization. The electromagnetic 
components vanish and the solution reverts to the original when a = 0. 

In terms of the notation developed in Chapter 16, the initial solution 
may be described either by the function Xo with uj = 0, or by ^ with 
<E> = 0. The new solution is then given by 

X = Xo(cos 2 a + sin 2 ae- c/ °Xo) , w = 0, (18.3) 

and the potential H can be obtained by integrating the equations 

H v = i sin a cos a (U ov + V ov ) e~ L 
H u = -i sin a cos a (U ou + V ou ) e 

Since H is purely imaginary, the new solution may also be described in 
terms of the potentials 

* = * (cos 2 a + sin 2 a e -2 ^- 1 ) 2 , $ = 0. (18.5) 

Hence the potential Z may be determined using (16.19). 

In his original paper, Panov (1978) applied this transformation to 
the Khan-Penrose solution. In his second paper (1979a), he applied it to 
the general solution given by Szekeres (1972) as described in Section 14.2. 
In these solutions there is always a space-like curvature singularity on the 
surface given by e~ u = f + g = on which the opposing waves mutually 
focus each other. 

The technique described in this section will be referred to again in 
the following section. It will be shown there how other solutions that 
do not contain strong curvature singularities, such as the Bell-Szekeres 
solution, may be included in the notation of this section. 

18.2 An alternative approach 

The class of solutions in which the metric can be diagonalized is charac- 
terized by the fact that the Ernst potentials Z and E are real and that 
uj = 0. In terms of the functions of the Szekeres line element (6.20), 
the equivalent condition is that W = 0, and hence from (6.22e), that 
$ $ 2 = $ 2 5>o- In the general case, the components $ , $2 and H or r/ 
must remain complex. 



(18.4) 
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In this case the main field equations (16.6) putting x = e ~ V > or 
(6.22 d, e) expressed in terms of the potential H, together with Maxwell's 
equations (16.7), can be written in the form 

2V UV - U U V V - U V V U + 2e u - v (H u H v + H U H V ) = 

H U H V - H U H V = (18.6) 
2H UV — V U H V — V V H U = 

In the particular case when the polarization vectors of the approach- 
ing gravitational waves and electromagnetic waves are all aligned, it is 
possible to use a duality transformation to ensure that $o, $2 and thus 
H or rj are all real. 

In order to generate exact solutions of this type, it is appropriate to 
consider the technique described at the end of Section 17.4. Accordingly, 
we consider an initial solution of the vacuum field equations describing 
the collision of plane gravitational waves. This may be described in terms 
of the initial Ernst potential E Q which is now real. The potentials for 
a new diagonal solution describing the collision of plane electromagnetic 
waves, usually coupled to gravitational waves, are then given by (17.32) or 
(17.34), where a, 6, c and d are all real and satisfy (17.33). Such solutions 
clearly revert to the original vacuum solution when a = 1, b = 0, c = 
and d = 0. 

It is of interest initially to consider the restricted class of solutions 
that can be obtained from (17.32) when 6 = and d = 0. In this case 
c = y/l — a 2 and the transformation reduces to (17.24). Then, reverting 
to the alternative parameter a given by (17.25), the new solution is given 
partly as in (17.27) by 

sin 2 a + cos 2 a Z Q sin a cos a (Z Q — 1) 

= ~ o = o 2 

cos 2 a + sin a Z Q cos 2 a + sin a Z Q (18.7) 

* = Z - H 2 = Z Q (cos 2 a + sin 2 a Z )~ 2 

The metric functions V or x are thus given by 

e~ v = x = e~ v ° (^cos 2 a + sin 2 ae- (c/o - 1/o) ) 2 . (18.8) 

This expression can be seen to be indentical to that of (18.1) apart from 
an insignificant change in the sign of V. It can therefore be seen that the 
technique described by (17.24) in the diagonal case is identical to that 
of Panov (1978, 1979a) that has already been described in the previous 
section. The complete solution in this case is given exactly as in (18.1) 
and (18.2) but with a change of sign for V and V Q . 
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It has been noted in the previous section that solutions obtained 
using the general technique (17.32) automatically satisfy the boundary 
conditions for colliding electromagnetic and gravitational waves if those 
conditions are satisfied for the initial vacuum solution. In the particular 
case of Panov, as described in Section 18.1, this has already been noted. 

It is important to notice that we are here using the Ernst potential 
rather than the metric functions. Thus using (11.3), (11.11) and (12.36), 
the functions associated with the initial solution are related by 

^ = 4 = l „ = iEME? = ( ^. (18 . 9) 

Using these identities, the main field equations may be seen to reduce 
to the single real linear equation 

((1 - t 2 )V , t ), t - ((1 - z 2 )V , z ), z = (18.10) 

which may be seen to be identical to (10.13), which is the main equation 
for vacuum solutions with aligned polarization. A general solution of this 
equation will contain the terms (10.16), (10.20) and (10.22), and can be 
written in the form 

V = Y,[anPn(t)Pn(z) + q n Q n (t)Pn(z) + PnPn(t)Qn(z) 

n 

+ b n Q n (t)Q n (z)} - |clog(l - t 2 ) - §clog(l - z 2 ) (18.11) 
+ ^ di cosh -1 | ; | 

where P n (x) and Q n (x) are Legendre functions of the first and second 
kinds respectively, and a n , p n , q n , b n Ci and di are sets of arbitrary con- 
stants. 

In his first paper, Panov (1978) derived a new solution using the 
above technique with 6 = and d = and with the initial solution taken 
to be the Khan-Penrose solution given as in (10.18) by 

V = -2Q (t)P (z). (18.12) 

In his second paper (1979a), he applied the same technique to the general 
solution given by Szekeres (1972). 

It may also be observed using (16.31) that the Bell-Szekeres solution 
can be obtained with this method by setting a = 0, 6 = 0, c = 1 and 
d = and using the initial solution of (18.8) given by 

V = -2P (t)Qo(z) - \ log(l - t 2 ) - \ log(l - z 2 ) 
= -2P (t)Q (z) + U . 
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However, as noted in (16.32,33), the same solution is also obtained if the 
term 2Po(t)Qo(z) is replaced by 2Qo(t)Po(z). It can then be seen from 
(10.35) that the initial solution, in this case, is the degenerate Ferrari- 
Ibahez solution described in Section 10.5. 

In their substantial paper that has already been frequently quoted, 
Chandrasekhar and Xanthopoulos (1987a) have also described a general 
class of diagonal solutions. This has also been obtained using the above 
approach but with the constraint that a = 0, 6 = 0, c = 1 and d = so 
that Z = 1. In this case, it is convenient to introduce a new function F 
such that 

F=\{V -U ). (18.14) 
With this, the new solution can be characterized by the potentials 

Z=l, # = tanhF (18.15) 

and the metric and field components can conveniently be expressed in the 
form 

U = U Q , 
W = 0, 

$° = -sechFF, 

Chandrasekhar and Xanthopoulos (1987a) have also obtained the 
interesting result that if, using the transformation (18.16), the initial so- 
lution is chosen such that 

oo 

F = a n P n (t)P n (z) - Q (t)Po(z) (18.17) 

n=0 

then, in the new solution, the surface on which t = 1 does not correspond 
to a curvature singularity as in the Bell-Szekeres solution and the solution 
of Section 17.3. 

Using the same definition (18.14), it is also possible to show that the 
more general solution obtained from (17.32) can be written in the form 

U = U , e~ v = .- 1 , e- u °((l -5)coshF-asinhF) 2 , 

(a 2 + c z ) v ' 

W = 0, e~ M = e 2F ((l-6)coshF-asinhF)V M °, 

$D _ (c + ad-bc)F v ^ D _ (c + ad-bc)F v 

" ~(1 -6)coshF-asinhF' 2 ~ (1 - b) coshF - asinhF' 

(18.18) 



e -v = e -^ cos h 2 F, 

e - M = e 2F cosh 2 Fe- M °, (18.16) 
$2 = sech FF U . 
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It may be noted that a very general solution may be obtained by 
taking the initial solution in the form (18.11), and using this generation 
technique given by (18.18) with F given by (18.14). It would appear 
that such general solutions would normally contain a scalar curvature 
singularity on the surface t = 1. However, although this class of solutions 
seems to be very general, it does not necessarily contain the complete class 
of all possible diagonal solutions, even for the restricted case in which $o 5 
$2 and thus H or r\ are all real. 

18.3 Electromagnetic Gowdy cosmologies 

It has already been pointed out in Section 10.7 that the Gowdy cosmolo- 
gies, which contain interacting plane waves, can under certain conditions 
also be made to satisfy the necessary boundary conditions appropriate 
for their interpretation in terms of a collision of initial plane waves. The 
extension to the electromagnetic Gowdy cosmologies has been carried out 
by Charach (1979). A reinterpretation of these solutions as the interaction 
following the collision of plane electromagnetic waves is now presented in 
the above notation. 

It is first appropriate to re-introduce the alternative coordinate sys- 
tem (10.58) defined by 

t = f + g, z = f-g (18.19) 

where t is a past pointing, time-like coordinate which decreases towards 
the singularity when t = in region IV. With this, the line element for the 
diagonal metric is given by (10.59), and the main field equations (18.6) 
take the form 

V +lv -V" + le~ v (HH - H'H') = (18.20a) 

t t 

HH' - H'H = (18.206) 
H — H" — VH + V'H' = (18.20c) 

where the dot and prime denote derivatives with respect to t and z re- 
spectively. 

To express these equations in the form of Carach (1979), it is conve- 
nient first to put 

F = 2^-logt, or ^=\{V-U). (18.21) 

It is then possible to put 



H = -p + le 2 V, H' = -p' + le 2 ^& (18.22) 

t t 
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where p and a are the two non-zero components of the four-potential for 
the electromagnetic field. With them the non-zero components of the 
electromagnetic field tensor are given by 

F±4 = —F41 = a, -F24 = — F42 = o~' , 

(18.23) 

F13 = — F31 = p, F23 = — F32 = p ■ 

With these substitutions, the main field equations together with 
Maxwell's equations take the form 

i> + \i> - r = -e-*V - p' 2 ) + l 2 e^(a* - a' 2 ) 
11 

pa — p a 

..1. „ . /: (!8-24) 

P+ r -p 

a -la- a" 



It may be noticed that the imaginary part of (18.20c) is automatically 
satisfied, and that the fourth equation in (18.24) enters as the integrability 
condition for the definitions (18.22). 

Charach has presented a family of exact solutions of these equations 
for which and p are functionally related by the equation 

V> = |logP(p), and (7 = 0. (18.25) 

In this case, H = — p, and the equations (18.24) become 








2(V>p- 


- W) 


2(^d 


- VV) 



(18.26) 



It may be assumed that p 7^ ±p' as this does not lead to solutions with 
interacting waves. It follows then that P can be written in the form 



P = ci + c 2 p - p 2 (18.27) 

where c\ and c 2 are arbitrary constants. The non-linearity in the resulting 
equation can then be removed by the substitution 

p = p - atanhaF (18.28) 



18.3 Electromagnetic Gowdy cosmologies 175 



provided 

Cl = a 2 - pi c 2 = 2 Po (18.29) 
and F is a solution of the linear equation 

F + If - F" = 
t 

It may immediately be noticed that (18.30) is identical to the vacuum 
equation (10.60) as may also be deduced from (18.20a). This equation 
may also be written in terms of different coordinates in the form of the 
vacuum equations (10.2) or (9.3) as 

2F UV - U U F V - U V F U = 0. (18.31) 

The resulting metric function is then given by e~ v = e~ u a~ 2 cosh 2 aF. 

The method described above can clearly be used to generate electro- 
magnetic solutions from any initial vacuum solution. It may be noticed, 
however, that the constant p can be removed from the electromagnetic 
potential, and that it is possible to set a = 1 without loss of generality. 
In this case the solutions are characterized by the potentials 

Z = l, H = tanhF (18.32) 

which can be seen to be identical to (18.15). This solution-generating 
technique is thus identical to that of Chandrasekhar and Xanthopoulos 
(1987a) as described at the end of the previous section. In terms of the 
function F, the metric and field components are given by (18.16). 
Charach (1979) has presented the class of solutions in which 

F = - alogt + ^{Au cos[u(z + a u )]J (ujt)} 

+ Y,i B " cos M 5 + 0u,)]Y o (wi)} 

where Jo(nt) and Yo(nt) are Bessel functions of the first and second kinds 
of zero order, and A n , B n , a n and j3 n are sets of arbitrary constants. This 
has been obtained by separating the variables in (18.30), and considering 
only solutions that are periodic in z. 

As pointed out by Feinstein and Ibanez (1989) and described here 
in Section 10.7, these solutions can easily be adapted to the boundary 
conditions for colliding plane waves if F also includes the terms 

F = ^^cosh" 1 ■ (18-34) 



(18.30) 



(18.33) 
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It may also be noted that, since the equation (18.30) or (18.31) for F 
is the by now familiar Euler-Poisson-Darboux equation, further explicit 
solutions can be obtained using any combination of the terms described 
previously in Section 10.8. 

18.4 Other methods 

Further diagonal solutions which describe colliding plane electromagnetic 
waves can easily be generated from known solutions, and also from known 
colliding plane gravitational wave solutions, using the invariance trans- 
formations (17.35a, d,e) with real constants. These transformations are 
generalizations of those such as (18.7) that have been considered above. 
Other methods are also available. 

Within the context of stationary axisymmetric solutions, it is well 
known that a Bonnor transformation (Bonnor, 1961) can be used to relate 
a class of diagonal electromagnetic solutions to a class of non-diagonal 
vacuum solutions. The precise result may be stated in the form that if Z 
is a complex solution of the vacuum Ernst equation (11.8), then a solution 
of the Ernst equations (16.22) for an electromagnetic field is given by 

Z = Z Z , H=\{Z -Z ) (18.35) 

where it may be noted that here H is purely imaginary. The inverse 
transformation also holds. 

It was first suggested by Halilsoy (1983) that this transformation 
could be used to generate colliding plane wave solutions. However, there 
is considerable difficulty with the application of the boundary conditions. 
If Z satisfies the junction conditions (11.26) for colliding non-aligned 
gravitational plane waves, then (Z,H) given by (18.35) do not necessar- 
ily satisfy the junction conditions (16.28) for colliding aligned electro- 
magnetic waves. All is not lost, however, since it is possible to start with 
solutions of (11.8) or (16.22) which can not be interpreted as colliding 
plane waves, and then to use the transformation (18.35) or its inverse to 
generate solutions which may be so interpreted for certain values of the 
free parameters. 

In view of the large number of solutions of Ernst's equation that are 
now known, such an approach to the generation of colliding plane wave 
solutions is feasible. However, it has not yet been used effectively. 



